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Abstract 

We develop an approach for performing scaling analysis of A-step Random Walks (RWs). The 
mean square end-to-end distance, {R%), is written in terms of inner persistence lengths (IPLs), 
which we define by the ensemble averages of dot products between the walker’s position and 
displacement vectors, at the j-th step. For RW models statistically invariant under orthogonal 
transformations, we analytically introduce a relation between {Pj^) and the persistence length, 
Aat, which is defined as the mean end-to-end vector projection in the first step direction. For 
Self-Avoiding Walks (SAWs) on 2D and 3D lattices we introduce a series expansion for Xjy, and 
by Monte Carlo simulations we find that Aoo is equal to a constant; the scaling corrections for Xjsf 
can be second and higher order corrections to scaling for (R/f). Building SAWs with typically one 
hundred steps, we estimate the exponents uq and Ai from the IPL behavior as function of j. The 
obtained results are in excellent agreement with those in the literature. This shows that only an 
ensemble of paths with the same length is sufficient for determining the scaling behavior of {R%), 
being that the whole information needed is contained in the inner part of the paths. 
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I. INTRODUCTION 


Random Walk (RW) models are ubiquitous in the literature with applications in several 
areas, such as Physics pQ, Biology [2] and Economy |3]. The simplest case is the walker 
displacement in a sequence of independent random steps, namely ordinary RW [1]. One 
may also obtain random paths on a geometrical space with distinct displacement schemes, 
leading to other RW models. A fundamental importance of these models lies in the fact 
that many real phenomena can be mapped or directly represented by paths traversed by 
walkers in some geometrical space, e.g., a single-strand DNA [3] and magnetic systems [3]. 
An example is the Self-A voiding Walk (SAW) dehned by a walker forming a random path 
that never intersects itself; standard SAWs are performed on regular lattices, where the 
walker steps to nearest-neighbor sites and does not visit a site more than once [7]. 

Because of non-overlapping paths, the SAW model plays a central role in Polymer Physics 
[8] by capturing the excluded volume effect in a dilute solution under good solvent condition 
or at high temperatures The SAW model is also well known in statistical physics context 
because of its equivalence with the n-vector model with n —>■ 0, as de Gennes hrst pointed 
out [To]. From this equivalence, with arguments of renormalization and held theories, one 
expects the following series expansion for the mean square end-to-end distance [mill]: 

{R%)m = aoAr2-o(i + ^ + ^ + ... 

+ + ( 1 ) 

where z/q is the leading exponent. The terms proportional to N~^ with i = 1, 2, • • •, are 
analytical corrections, and the terms proportional to with non-integer exponents 

Aj < Aj+i and j = 0,1, 2 ■ • •, are the non-analytical corrections to scaling. The leading and 
corrections to scaling exponents are universal. The indexed brackets {.)n refers to the N- 
step RW ensemble average, and from now on, unless strictly necessary, we omit the index N. 
Numerical estimates of exponents z/q and Ai are based on either exact counting techniques 
psun, or in Monte Carlo (MC) simulation methods [131 [131, through the sampling of {R%) 

[muH]. 

Obtaining such estimates for uq and Ai, especially for 3D SAW, is a challenge from 
several points of view. The exponential growth of the number of possible A-step paths 
cn ~ where /i is the connectivity constant and 7 > 1, imposes a limit to exact 
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counting. To the best of our knowledge, the maximum values obtained are N = 79 [19] 
and N = 36 [T3| for SAWs on 2D and 3D square lattices, respectively. Concerning Monte 
Carlo simulations, there exist an appeal to hnd z/q and Ai using very long paths. Obtaining 
high quality Monte Carlo data for such path lengths is an extremely difficult task for the 
SAW model. The variable length algorithms suffer from attrition problems, namely barriers 
that prevent paths to grow, while the hxed length algorithms suffer from the decreasing of 
acceptance rate to generate a new non-self-intersecting path, according to the increase of 
the (hxed) path length [20]. 

Numerical drawbacks also take place when one studies other conformational quantities. 
An example is the persistence length, Xn, dehned as the mean end-to-end vector projection 
in a hxed direction along the hrst step [211 [22], as A —)■ oo [23]. Dehning the end-to-end 
vector as Rn = where Uj is the walker displacement at the j-th step, the persistence 

length can be expressed by Aat = {Rf^-Ui)/\ui\ = Numerical results of A at, 

for 2D-SAWs, are controversial in the literature, and for 3D, are scarce [32] • For 2D-SAW, 
Grassberger [33] obtained the hrst estimate of Aat in the square lattice, by means of a power 
law Atv ~ , with 6 = 0.063(10). Since for 0 0, it is also well htted by Aat ~ hi(A), 
as suggested by Redner and Privmann [33] • They obtained both estimates by sampling 
the displacements projections along the hrst step direction, for all possible conhgurations of 
SAW paths with N < 24. This weak divergence has been questioned recently by Eisenberg 
and Baram [33], because their MC estimates of {ui -Uj) show that Xn converges to a constant 
when N —)■ oo. One could employ Aat in Monte Carlo [2H] and experimental characterization 
of certain polymers [291 ED], despite there exist some limitations of Aat measures such as 
divergence and edge ehects [31] . 

Rehned results about the scaling behavior of the aforementioned conformational quanti¬ 
ties to study universality are challenging, and have been the subject of discussion for many 
years [HIITS]. As usually one does not have exact results for the SAW model, there exists an 
appeal for simulations of large, sometimes very large, paths. Here, one proposes to answer 
two questions about a SAW: (i) What is the asymptotic limit of its persistence length? (ii) 
Is there some way to hnd out its scaling behavior employing relatively small chains? To 
answer these questions, we found an approach for performing scaling analysis of RWs, by 
focusing in the behavior of {R%)- 

The structure of the paper is as follows: In Sec. [IT] we present the analytical results by 
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defining the inner persistence length and their relation with {R%) and Xn, for RW models 


statistically invariant under orthogonal transformations. In Sec. we provide a series 
expansion for Xn and obtain the scaling behavior of 2D and 3D-SAW models with Monte 
Carlo simulations; we also obtain reliable estimates of the exponents z/q and Ai and discuss 
the contribution of Ajv to {R%) behavior. In Sec. 


IV 


we give concluding remarks. 


II. INNER PERSISTENCE LENGTH AND ANALYTICAL RESULTS 

We dehne the inner persistence length (IPL) for an A-step RW, by the average dot 
product: Xj = {Rj ■ Uj). To relate {R%) to Ij, and X^v to X^, we write the square distance 
at the j-th step for an A-step RW as: R^ = R^_i + ^Rj ■ Uj — v?-. Adding up we have 

Ej=i R] = Ei=i R]-i + Ej=i 2 Ri ■ Uj - Y!i=i Ifof, where Ro = 0 leads to Y!}=i R]-i = 
Thus, considering \uj\ = 1, we write the average {R\) = ‘^Yl^j=i'Rj ~ fo 
particular for k = N, the mean square end-to-end distance is 

N 

( 2 ) 

t = l 

Now, consider a generic class of RWs, where ensembles of A-step walks obey the following 
invariance property: the probability distributions, of each step fo, i = 1,2, ...,A, which 
compose a path, is invariant under orthogonal transformations. With this, we exclude walks 
like the tourist model [32], where the medium disorder ra breaks down such invariance 
symmetries. Particularly, one considers an ensemble of A-step RWs obeying the mentioned 
probabilistic symmetry, under a specihc orthogonal transformation T given by fo —>■ 
the prime denotes the displacement vectors in the transformed reference frame, and u[ = 
—UN-i+i, with i = 1,2,...A. Notice that u' G {ui,'U 2 ...'Uv}, where {.} represents the 
complete ensemble of paths. This symmetry operation can be achieved by a translation 
followed by inversion of all displacement vectors. In other words, one does invert each path 
and change the origin to the end of the walk. An immediate consequence for the complete 
ensemble of random paths is {fo} = {Wi}, with z = 1,2, ...A, which leads to {Rat} = {Rjv}- 
From the previous relations, it follows that {Rat • un} = {Rn ■ fo}, so the conhgurational 
average (Rat ■ mtv) = {Rn ■ fo) holds. This average, for A —)■ cx), is the persistence length 
A at. Therefore, the mean square end-to-end distance could be rewritten as 

(R^) = (R^_i) 2Xis[ — 1, (3) 
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and we have established a relation between {R%) and Xn- We observed Eq. nnmerically, 
prior to its proof, by exact calculations for N < 24. Some RW models that obey such a 
relation are the iV-step ensemble of ordinary RW and SAW paths. 


III. NUMERICAL RESULTS FOR THE SAW MODEL 

From now on, we numerically study Xj for SAWs using the non-reversed random walk 
(NRRW) algorithm to generate the ensemble of A-step non-overlapping paths. Because 
of the attrition problem, i.e., barriers or traps that prevent paths to achieve N steps, the 
NRRW is inefficient to generate good statistics for long SAWs, since the probability decays 
as ptv cc exp[— 7 A], where 0 < 7 < 1 is the attrition constant. However, the generated data 
with this algorithm are surprisingly good enough to validate our approach, showing that we 
choose the right corrections to scaling terms in the expansion of IPLs. 

Starting with {R%), we now analyze the persistence length. For the square lattice, uq = 
3/4 [3H] and a common belief is that Ai = 3/2 [30] • With these exponents values, from 
Eq. using only the hrst two leading exponents, we see that {R%) ~ AN^R + BN^R. 
The same reasoning leads to a similar result for cubic lattices, where vq ~ 0.587597(7) and 
Ai ~ 0.528(12) are widely accepted values [Hj. Both averages in Eq. {R\[) and {R%_i), 
are obtained considering the same A-step ensemble. In this sense, we follow our previous 
notation by omitting the bracket index. The difference {R%) — {R%_i) seems to be the 
discrete derivative of square end-to-end distance, which is not true for the SAW model. 
One should evaluate the derivative considering SAW ensembles of A and (A — l)-steps: 

— {Rn-i)n-i- According to Eq. the leading term of {R%)n derivative is 
with the hrst two corrections proportional to and respectively. From the 

persistence length plots in Fig. Atv clearly does not diverge as the leading term of {R%) 
derivative, instead it seems to converge to a constant as A goes to inhnity [SH]- Thus, we 
introduce the following series expansion: 

= ao + + asA^""" + • • • (4) 

where the exponents Wi > 0, i = 1,2,3---, are linear combinations of uq with analytical 
and non-analytical corrections to scaling exponents. As for example, from the persistence 
length data htting with Eq. (see Fig. [^, we hnd that Wi = 2 uq — 2 and W2 = 2uq — Ai — 1 
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TABLE I. Coefficients and exponents for fitting, with Eq.|^ the Atv data obtained from simulations 
for 2D and 3D square lattices. The wi = 0.34(5) value is an effective exponent, thus depending on 
the coefficients ai and exponents Wi of Eq. [40] . 


d 

oo 

ai 

02 

Wl 

W2 

2 


2.525(4) 

-2.32(3) 

0.81(3) 

0.5 

1 


1 

2.664(3) 

-1.714(9) 

- 

0.34 

- 

3 


1.422(1) 

-0.39(6) 

-0.022(5) 

0.8248 

0.34 


Fitting with equation Aw ~ oq + uiA 0-34(5) from Ref. [35] . 


are the best choices. The ctj and Wi values are shown in Tab. [Tj An immediate consequence 
of such hndings along with Eq. is that Aw could contribute only with second and higher- 
order of analytic and non-analytic corrections for {R%). Our estimate of Aw, for square 
lattices, is compatible with the one of Eisenberg and Baram [35]. Through their estimate 
of the step-step correlation scaling: {ui ■ Uj) = (^ij)w ~ and the dehnition 

Aw = j)w, we obtained Aw ~ cto ~ 1-7N with which we htted the persistence 

length data, but leaving ai free, as shown in the inset of Fig. [^a). 




FIG. 1. SAW persistence length for (a) square and (b) cubic lattices. In both lattices, Aw 
converges to a constant. The inset plot (a) shows Aw fitted by the function ao + of 

Ref. [35], where Aqo = 2.664(3) is compatible with our estimate Aoo = 2.525(4) from Eq. For 
the cubic lattice, Aqo ~ \/2 is compatible with the one of Ref. [32]. The inset plot (b) depicts the 
random pattern of the residual plot for Aw when fitted by Eq. [^ 

Now, consider Ij, for 1 < j < Al. According to the collapsed log x log plots of Fig. it is 
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FIG. 2. IPL data collapse in log-log scale for (a) square and (b) cubic lattices with {A N = 30), 
(♦ = 40), {MN = 50) and {• N = 60). The {Rj ■ Uj)N behaves as linear increasing function up to 

~ jmax, with a slope « {2i/q — 1). In both lattices jmax oc N, with constant of proportionality close 
to each other (~ 0.7). For j > jmax the scalar products contribute to residual terms of corrections 
to scaling of {R%)- 

notable that Ij looks like a straight line up to near the point where it reaches its maximum 
value, at the jmax step, with a positive slope ~ 2z/o — 1. From Eqs. and and Fig. 
assuming Zj scales as jg reasonable, at least for j < jmax- Such proportionality leads us 

to look for reliable estimates of z/q, and corrections to scaling exponents, for SAW ensembles 
with N not too large. To accomplish this aim, diminishing the influence of the A^-step 
ensemble on estimates of scaling exponents is necessary. In other words, it is necessary to 
hnd a cutoff step j = jc{N), at which Zj begins to be noticeably influenced by the A^-step 
SAW ensemble. Surely, we can neglect steps above jmax- To seek the jc{N) step, we use the 
difference between the IPLs of two A^-step ensembles, one that contains A^i, and the other 
N 2 steps, 

ARj{Ni, N 2 ) = {Rj ■ Uj)N2 - {Rj ■ Uj)Ni, ( 5 ) 

where A '"2 > A^i- According to Fig. |^a), the IPL has approximately the same behavior for 
the two path lengths, up to the middle of the shortest path, jc(A'i) ~ A'i/2, for square 
lattices. Similarly, for cubic lattices, it has the same behavior, up to a third of the shortest 
path jc{Ni) ~ Ni/3 [see Fig. |^b)]. Therefore, using j < jc{N), with jc{N) = N/2 and 
jc{,N) = A^/3 for 2D and 3D lattices, respectively, it is suitable to estimate the scaling 
exponents through Zj. 

Additional information to do scaling analysis with Zj comes from the expansion of {R%) 
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FIG. 3. IPL differences {ARj{Ni, N2)) for SAWs: (a) for square and (b) cubic lattices, with 
Ni = 40 and N2 = 60, and A^i = 60 and N2 = 90, respectively. According to ARj{Ni, N2) 
depicted here [see Eqj5], the Ij starts to be influenced for jc > Ni /2 and jc > Ni/ 3 , for 2D and 3D 
square lattices, respectively. Inset plots show IPL non-weighted fit using Eq. within a confidence 
interval of 95%. The square lattice data includes N ranging from 90 to 120 with increment of 10. 
The fitting parameters obtained are tp = 0.4979(21) and /3i = 0.6630(50). The cubic lattice data 
includes N ranging from 150 to 195 with increment of 15. The fitting parameters obtained are 
ip = 0.1752(14), Ai = 0.522(52), and /3i = 0.7581(56). 

in powers of N . We have found no evidence of the linear term in the expansion of {R%) on 
square or cubic lattices. The nonexistence of the linear term is also reported in Refs. mills]. 
From Eq. the only way to disappear with the linear term in the expansion of {R%) is if 
the summation of Xj cancels it out. This Ending, with Eq. leads us to write 

X,- = /3o + /5i(i - tY [1 + P 2 U - r)-^^ + (6) 

for i < jc{N), where r is a smoothing constant [13|- We set (3o = 1/2 just to cancel the 
linear term. Also, we did another ansatz: (32 = — (2i/o — 1) and r = 0.5. This was inspired 
by results considering only the first non-analytical correction to scaling term, and leaving 
only the parameters /di, (32 and r free, which lead us to And (32 ~ — (2i/o — 1) for the 3D case. 
Notice that, in general (3i = (3i {N) and (32 = (32 (N); however, for N not too large, order of 
hundreds for 2D and 3D cases, these parameters converged to constants, for j < jc{N). 

The IPL data, containing several A^-step ensembles, fitted by Eq. is depicted in the 
inset plots of Fig. For both, the 2D and 3D square lattices, the leading and sub-leading 
exponents are in excellent agreement with the believed results. For the square lattice, we 
found z/q = 0.7489(21), and the non-analytical hrst exponent results in Ai = 3/2; because it 
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does not appear in Eq. showing that there exists a constant in the expansion of {R%)- This 
is confirmed through the expansion of Xn] the predicted results are ip = 0.5 and Ai = 3/2. 
For cubic lattices we found z/q = 0.58757(140), and Ai = 0.522(52), while the best predicted 
results are z/q = 0.587597(7) and Ai = 0.528(12) [T7] . Using several A-step ensembles 
seeks to reduce the error on exponent estimates; however, they may carry some small biased 
errors. To check this, for 2D-SAW, we used N = 120 steps obtaining uq = 0.7500(63), and 
for 3D-SAW we used N = 198 steps giving z/q = 0.58758(450) and Ai = 0.52(17). However, 
the errors we get are not as small as those from literature for the 3D case mi- We can 
improve these results, by taking into account the advantage of the statistical invariance, and 
calculating the IPL starting from the end of the generated chains, thus doubling the sample. 
In fact, it is out of the scope of this paper to find high precision values for the exponents, but 
to validate and evaluate the benefits of our approach. Moreover, the whole potential of the 
method to do the scaling analysis of RWs has not been fully exploited. We expect that the 
corrections to scaling exponents are easily accessible from the study of the monotonically 
decreasing Zj terms of {R%), which will readily be tackled. 


IV. CONCLUDING REMARKS 

In summary, we have proposed an approach to address the scaling of RW conformational 
quantities, where the mean square end-to-end distance is proportional to the summation 
of the inner persistence length, Zj = {Rj ■ Uj). For RW models, where paths obtained by 
orthogonal transformations occur with the same probability, we obtained a novel relation 
between the mean square end-to-end distance and persistence length. Despite the numerical 
limitations to do scaling analysis, we introduce a series for the persistence length Aat and show 
that it converges to a constant, ao, apart corrections to scaling terms. We also developed 
a method to calculate the scaling exponents from Zj with a path cutoff that diminishes the 
A-step ensemble influence. Thus, the method is efficient to obtain the scaling behavior of 
SAW. 

We conclude that only an ensemble of paths with the same length is sufficient for per¬ 
forming scaling analysis, being that the whole information needed are contained in the inner 
part of the paths. The scaling method discussed in this paper can be important for studying 
universality, criticality, and conformational properties of systems mapped on RW models. 
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such as polymers, biopolymers, and magnetic systems. 
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